Kondo effect in Complex Quantum Dots in the presence of 
an oscillating and fluctuating gate signal 
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We show how the charge input signal applied to the gate electrode in a double and triple quantum 
dot may be converted to a pulse in the Kondo cotunneling current being a spin response of a nano- 
device under a strong Coulomb blockade. The stochastic component of the input signal results 
in the infrared cutoff of Kondo transmission. The stochastization of the orbital component of the 
Kondo effect in triple quantum dots results in a noise-induced SU(4) — > 5(7(2) quantum transition. 
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I. INTRODUCTION 

Current interest in charge-spin conversion effects is 
spurred by challenging prospects of spintronics. Most 
of the mechanisms of such a conversion are related to the 
interconnection between electrical and spin current due 
to spin-orbit interaction^ which results in spin accumu- 
lation near the sample edges. Such an accumulation in 
three- and two-dimensional electron gas in elemental and 
III-V semiconductors may result in a spin-Hall effect 2 
and positive magnctoresistance£ It was argued also that 
the Rashba-type spin-orbit interaction in a quantum dot 
assists pure spin current by modulation of the voltages 
applied to the leads in a three-terminal device^ Spin 
Coulomb drag effects should be mentioned in this con- 
text, which may result in spin polarization of charge cur- 
rent due to intrinsic friction between electrons with dif- 
ferent spin projections induced by Coulomb scatteringi^ 

In all these propositions the possibilities of conversion 
of charge current into spin current were discussed. It is 
possible also to try to use the external electric field for 
the generation of spin current or another spin response. 
One such idea was formulated recently for light emitting 
diodes (LED) based on conjugated polymers,— where dis- 
sociation of excitons in a strong enough electric field may 
result in the accumulation of up and down spin densities 
near the two ends of the LED. 

In this paper we show that the charge input signal ap- 
plied to the gate electrode in double (DQD) and triple 
quantum dots (TQD) forming closed loops (equilateral 
triangle) may be conversed to a pulse in Kondo cotun- 
neling current, which is in fact the spin response of DQD 
under strong Coulomb blockade. The general idea of such 
a conversion was formulated in our previous paper— (here- 
after referred as I). Using the example of a T-shaped 
DQD occupied by two electrons we have shown that 
the charge-spin conversion is possible due to the spe- 
cific dynamical symmetry of a spin multiplet consisting 
of two singlets and a triplet. Since the states in this 
multiplet are constrained by Casimir operators of the 



group SO(5) characterizing the dynamical symmetry^ 
the time-dependent perturbation in the charge sector 
which affects only the singlet states results in a time- 
dependent potential acting on the triplet states. 

An important aspect of the problem is the interplay 
between the coherent (adiabatic) and stochastic com- 
ponents of the input signal. In paper I we concen- 
trated mainly on the coherent signal and discussed in 
details the conversion of charge noise into spin response 
in Kondo tunneling. This response may be interpreted 
as dephasing and decoherence of the Kondo screening. 
Here we develop the general scheme, where both coher- 
ent and stochastic response of Kondo tunneling to the 
time-dependent gate voltage are calculated. Besides, the 
decoherence due to the noise component of the input sig- 
nal is considered in the long relaxation limit contrary to 
I, where the white noise approximation was considered 
for the noise correlation function. 

Another mechanism of charge-to-spin conversion may 
be realized in an equilateral triangular TQD in an exter- 
nal magnetic field penetrating the triangle plane. This 
model was introduced in Ref. [13 with the purpose of 
studying the interference between Kondo tunneling and 
Aharonov-Bohm interference (see Ref. [ll| for an exper- 
imental realization). The actual dynamical symmetry 
of TQD in a contact with metallic electrodes is SU(4) 
because the system possesses two spin and two orbital 
degrees of freedom. We show in this paper that the time 
dependent gate potential affects only the orbital compo- 
nent of multiplet, but the Kondo tunneling is sensitive to 
this perturbation because both orbital and spin discrete 
degrees of freedom are involved in Kondo screening. 



II. COHERENT AND STOCHASTIC CHARGE 
SIGNAL 

The subject of our calculations is the study of the 
transformation of the charge input signal into a Kondo 
response of complex quantum dots in tunnel contact with 
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source and drain electrodes. We study the mechanism of 
activation of internal spin degrees of freedom by means of 
a time-dependent gate potential applied to passive elec- 
trode. Two examples of complex nano-devices will be 
considered. The first is an asymmetric double quantum 
dot (DQD) which contacts with metallic leads in a so 
called T-shape geometry (Fig. QJi). The second is a triple 
quantum dot (TQD) in the form of an equilateral trian- 
gle in a three-terminal configuration, where the bias is 
established between the dots 2,3 and the time-dependent 
gate potential is applied to the dot 1 ((Fig. [lb). 
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FIG. 1: (a) Double quantum dot in a T-shape two-terminal 
geometry (b) Triangular quantum dot in a three-terminal ge- 
ometry in a magnetic field perpendicular to the plane of the 
triangle. 



In both cases the electron tunneling is described by the 
Anderson Hamiltonian 

H = Hi, an d + Hdot + H tu n- (2.1) 

where three terms are related to band electrons in the 
leads, electrons in the complex dot and tunneling cou- 
pling between two subsystems, respectively. The lead 
Hamiltonian has the form 

Hb an d = y y^,£bk c bkcr C bka I (2-2) 

b 

Here b stands for the leads, k, a are the wave vector and 
spin projection, respectively. The leads are identical in 
our model, (sbk = £&) so it is convenient to rc-cxpand 
the lead electron states ct, in terms of irreducible basis 

OK(7 

states Cp ka of the corresponding point symmetry groups 
(mirror symmetry for DQD and triangular symmetry for 
TQD). 

The Hamiltonian of the dot is written as 

H dot £ //': • Y. H jf + F i W » ( 2 - 3 ) 

3 jj' 

The potential wells in a complex QD are enumerated by 
the index j. Here 

Hj = £ 3 n 3 + Q3»r 



describes electron states in the potential well j under the 
Coulomb blockade Qj , and rij are the occupation number 
operators. The second term 

stands for the interdot tunneling. By convention j = 1 is 
reserved for the well coupled to the gate. This coupling 
is described by the last term in (12.31) . 

#i(t) = \y g (p) + v g (t)]m , 

where the gate voltage contains both a static component 
V g (0) and a time-dependent perturbation v g (t) . 

The tunneling term in the Hamiltonian (|2.ip has the 
form 

H tun = wY J Y.^lkJ ] a + H.c.) 1 (2.5) 

The values of j and (3 are determined by the geometry of 
the complex QD (see below). 

We consider a general situation, where the gate poten- 
tial applied to a multi-valley complex QD contains both 
coherent and stochastic components 

Vg{t)=V g {t)+SVg(t). (2.6) 

The v g (t) is the coherent (deterministic) contribution and 
8v g (t) is the stochastic noise component which is defined 
by its moments 

HW = (2-7) 

8Vg{t)5Vg(t') = rff{t - t') 

The ovcrlinc stands for the ensemble average, the char- 
acteristic function f(t — t') will be specified below. 

The starting point of our investigation is the canonical 
transformation which converts the gate potential into a 
time-dependent operator involving one of the generators 
of the group characterizing the dynamical symmetry of 
a complex QD^^iii This transformation applied to the 
Schrocdingcr operator —ih(d/di) + Hdot give o 8 ' 12 i 13 

H dot = U^HdotUt - ihUr 1 -^-, (2.8) 

with 

Ui = exp[i0i(i)m] (2.9) 
and the phase <f>i(t) given by 

Mt) = \[ dt'v g (t'). (2.10) 

One may apply the Hausdorff expansion to the first term 
in (Ml) 

HaoM = ffffi + £ [SrASuH^]]...]^.!!) 

771 
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where Si = 4>i(t)ni, and H d J t includes all time- 
independent terms from (|2.3p 

It is expedient to introduce " even" and " odd" hopping 
operators 

T i? = Y} d i^ ±d \M (2-12) 

(7 

To the lowest orders in V g , the time dependent part of 
the dot the Hamiltonian acquires the form 

SH dot (t) =~VJ2 (*Mt)T^ + ^WT^\2.13) 
where 

1 r* 

Mt) = fil dt'v g (t') (2.14) 

Thus, we obtain the effective time-dependent dot 
Hamiltonian 

#dot = Hfl + SH coh (t) + SH stoch (t), (2.15) 

where the time-dependent perturbation contains a coher- 
ent component [the first term in (|2.13p ] and a stochastic 
one [the second term in (|2.13p ]. In many cases the coher- 
ent perturbation may be easily taken into account in an 
adiabatic approximation^i^ whereas the second term in 
(|2.15[) results in the stochastization of the quantum state 
of the complex QD. 

We are interested in the influence of such a time- 
dependent perturbation on the electron cotunncling 
through complex QDs in Kondo regime. To investi- 
gate this influence one should derive the effective spin 
Hamiltonian from (|2.1[) by means of a time-dependent 
Schrieffer- Wolff (SW) transformation^*^ taking into ac- 
count the perturbation SHdot(t) in (|2. 15[) . The adiabatic 
component of this perturbation results in temporal oscil- 
lations of Kondo transparency and enhances the tunnel 
conductance on average^ whereas the stochastic compo- 
nent of the gate potential is detrimental for Kondo tun- 
neling. It results in the loss of coherence of the Kondo 
singlet state and in the smearing of zero bias anomaly 
in tunnel conductance. In accordance with the general 
approach to decoherence phenomena^ one should dis- 
criminate between the decoherence of the ground state 
of a quantum-mechanical ensemble and its manifestation 
at a finite energy /temperature. In the latter case one 
should argue in terms of dephasing due to elastic and 
inelastic scattering. Both processes are relevant in our 
system (see I). 

In the two next sections we will show how the coher- 
ent part of 8H{t) results in the conversion of the charge 
input signal to a Kondo response, and why its stochastic 
component brings an end to this process. 



III. DOUBLE QUANTUM DOT IN T-SHAPE 
GEOMETRY. EVEN ELECTRON OCCUPATION 

We start with the T-shaped DQD (Fig. Q}t) and con- 
sider the case of even occupation Af = 2 with one elec- 
tron per potential well. In this two-terminal geometry 
the irreducible set Cp ka consists of two combinations of 
source (s) and drain (d) leads. Only even standing wave 
Ceka = (c s ka + c d ka) / V% enters -fftun (|2.5p . and we omit 
in what follows the index /3 = e. The intradot indices 
jf have two values 1,2. We consider a DQD with two 
equivalent wells (Qj = Q) and assume that the time- 
independent component V g (0) of the gate voltage modi- 
fies the single-electron spectrum in such a way that the 
charge transfer excitation due to hopping Hf 2 (|2.4p is a 
relatively soft excitation with the energy A12 <§C Q (see 
Fig. Hi.) 
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FIG. 2: Energy levels of an isolated DQD. (a) Single elec- 
tron levels. The interdot tunneling V is shown by a dashed 
arrow; (b) Two-electron levels Ea; (c) Evolution of Ea with 
the growing of the scaling parameter £ = \nDo/D (see text 
for further details). 

Under these conditions the low-lying part of the energy 
spectrum of DQD in the charge sector J\f = 2 consists of 
two singlet states and one triplet stated 

E s = £1 + £2 - 2aV, 

E T = £l +£ 2 , (3.1) 

E E = 2e 2 + Q2 + 2aV. 

Here a = V/Aes is the effective indirect exchange pa- 
rameter, which favors antiparallel orientation of electron 
spins in two valleys of the DQD. The ground state of the 
isolated DQD is the spin singlet Es- 

In spite of this, Kondo tunneling through DQD in T- 
shape geometry is possible under certain conditions, be- 
cause the contact with lead electrons renormalizes the 
effective exchange in such a way that the singlet /triplet 
level crossing is possible. This result was obtained in Ref. 
[Til by means of a renormalization group (RG) technique^ 
with scaling parameter £ = ItiDq/D, where Dq and D is 
the initial and current energy scales for electrons in the 
leads, respectively. 

Having in mind this structure of low-lying states, it 
is convenient to represent Hdot m terms of Hubbard op- 
erators X AA = I A) (A' I, where |A) are the eigenvectors 
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corresponding to the eigenstates Ea (|3.1[) : 

A 

K = Tv 1 S 1 E, i/ = ±l,0 

It is important for further calculations that the system of 
operators describing transitions between the levels of any 
multiplet consisting of two singlets and one triplet forms 
a closed set of generators of the 50(5) algebra. Ten gen- 
erators forming this algebra are packed into three vectors 
S, P, M and one scalar A. Vectors describe transitions 
within the triplet Tv, between the triplet and the singlet 
S and between the triplet and the singlet E, respectively. 
The scalar A stands for transitions between the singlets 
5 and E. All these operators may be expressed via Hub- 
bard operators X AA (see Refs. I^lilliol for further details). 
One may rewrite the Hamiltonian H dot in terms of these 
generators 

ffdot = \ {E T S 2 + E S P 2 + E E M 2 ) + Q(JV - 2) 2 (3.3) 

Besides, the hopping operator (E32) ma y be repre- 
sented as 

= 1AV2, (3.4) 

Finally, the Casimir operator C for the 50(5) group is 

C = S 2 + P 2 +M 2 +A 2 = 4 . (3.5) 

Usi ng Eqs. (pT3]l . (|3T4|) as input data i n (pTTTj) and 
(|2.13p . we obtain the effective Hamiltonian (|2.15p in the 
following form (see I) 

H dot (t) = i (S T S 2 + E s (t)P 2 + E E (t)M 2 ) - fi(C - 4) 

(3.6) 

with 

Es(t) = E s — 6 coh (t) - <5 st och,s(0 

E E (t) = E E + S coh {t)+S stochiE (t). (3.7) 

Here coherent and stochastic corrections enter in the 
form of time-dependent " energy levels" : 

6 coh (t) = (2V 2 1 & ES )4>\{t) , 
<Wh,s(<) = (V 2 /4A ST )^(fj , (3.8) 
SstochMt) = (V 2 /4A ET )tf(t) , 

Aaa' = \E\ — E\r\. As was pointed out in Ref. H, a 
charge perturbation cannot directly affect spin degrees 
of freedom, and the triplet level Et remains time inde- 
pendent. 

It is seen from (|3.7p . that the coherent component re- 
sults in the time-dependent shift of energy levels. It may 
be treated as an adiabatic correction provided the time- 
dependent perturbation is weak enough, 5 co h(t) <C A E t- 
Below we adopt this adiabatic approximation. 



Unlike the coherent renormalization, stochastic com- 
ponents are not "in phase", i.e., S at0 ch,s(t) ^ <Wch,E(*)- 
This inequality makes the constraint imposed on spin dy- 
namics by the Casimir operator C (|3 . 5[) "fragile" for the 
50(5) group. Another source of stochasticity is the last 
term in Eq. (|2.8[) . In lowest order in the stochastic corre- 
lation functions (|2.7p its contribution is ~ Kn\<j)i(t)<f)\(t). 
This contribution may be converted into additional 4-th 
order corrections to Eqs. ()3.7|) . and we neglect it in the 
following calculations. 

A time-dependent SW transformation of the Hamil- 
tonian (|2.ip is performed by means of canonical trans- 
formation H = UiHU^ 1 . The phase T in the matrix 
U-2 = expiT is given by 

T W = E H{t)X Sr *c ka + v E {t)X E ^c ka - H.c] , 

ka 

(3.9) 

The coefficients v E (t) are fixed through the con- 

dition 

dT 

H tun + [T, (H dot + H band )] = itl—. (3.10) 

The solution of this equation is described in I, and the 
resulting cotunncling Hamiltonian has the form 

H cotun (t) = J T S-s + J s (t)P-s + J E (t)M- 43.ll) 

Again, pure spin scattering is not affected by charge per- 
turbation, but time-dependent spin-flip transitions in the 
leads described by the two last terms in p. lip arise due 
to the fact that the dynamical symmetry of the dot spin 
multiplet is activated by the time-dependent gate poten- 
tial. Like in Eq. (|3.6p . the time-dependent coupling pa- 
rameters in the SW Hamiltonian contain both coherent 
and stochastic components, 

J A (t) = 4 + J A d (t) + JL ch (t) (3.12) 

(A = S,E). The time-dependent corrections to Jq are 
calculated in Appendix A. Now we are well prepared to 
study the contribution of adiabatic and stochastic cor- 
rections to Kondo tunneling. 



A. Coherent input signal. 

We study here the transformation of the monochro- 
matic gate potential 

Vg(t) = VgCOSClt. (3.13) 

into a coherent (adiabatic) Kondo response under the 
condition fl <C A$t and <5 C oh(£) C A^. Then the phase 
4>\ (t) has a simple form 

M*) = ^(sinfli-l) (3.14) 
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We calculate the coherent (adiabatic) response at low 
enough energy and temperature, where the last term in 
p. lip may be omitted. Then we remain in the reduced 
{S, T} part of Hubert states, i.e. assume that the Kondo 
temperature Tk ~ Ast is valid. The dynamical symme- 
try of the reduced adiabatic (ad) effective Hamiltonian 



Hsw 



= H, 



(ad) 



dot 



H, 



(ad) 



(3.15) 



H, 



(ad) 



2 (^tS 
J T S • s 



2 f£ s (t)p 2 ) + g(^- 

[J S + jf d (t)]P-s 



2f 



is 50(4). 

The adiabatic part of the time-dependent Hamiltonian 
may be incorporated in a Haldane-type RG theor y 8 ' 12 ! 17 
As a result, the levels Ea acquire self energies scaled with 
the parameter £ = \u(Dq/D), namely Ma = o;a£ so that 



Ej 1 
Es(t) 



E T - a T £, 



(3.16) 



The self energy Ms depends paramctrically on time. The 
coefficients a s T were calculated in Ref. Il4 , and it was 
shown there that the inequality a T > a s is always valid 
due to existence of excited singlet level Ee- Due to this 
inequality, the S/T level crossing may occur at this stage 
of renormalization (see Fig. [2fc), so that the ground state 
of the system is triplet to the right of the crossing point. 
The time dependent factor as(t) describes parametri- 
cally slow variations of the scaling trajectory due to adi- 
abatic (coherent) corrections given by Eqs. Q3.7[) and 



The expression for J^ d (t) derived in Appendix A reads 



explicitly (see Eq. [A3 



J s 



ad 



V2W 2 ■ 
€ 2 -M s ' 



l(*) 



(3.17) 



(one may neglect the time-dependence of the denomina- 
tor). 

It is known from Kondo theory for quantum dots with 
SO (4) symmetry ; 14 i 20 i 21 that the Kondo temperature is 
a very sharp function of Ast with a maximum Tko a t 
Ast = (see Fig. [3j left column). It has an intermediate 
asymptotic behaviour for positive Ast where the ground 
state of DQD is a triplet 



Tko 



T 



KO 



A ST (t) 



(3.18) 



valid in an intermediate asymptotic regime for positive 
Ast at Tko/ Ast ^ 1- Here rj < 1 is a universal con- 
stant. This sharp dependence is a key to the charge- 
spin transformation mechanism, which is especially ef- 
fective in the vicinity of the triplet/singlet critical point 
Ast = 0. 

Next, we estimate the influence of the trembling signal 
on the tunnel conductance G(T, t) at given T > Tk in 
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FIG. 3: (Color online) Left column: Tk as a function of Ast- 
Middle column: time dependent 7V(t) corresponding to the 
evolution of Tk(Ast) in the left column. The intervals over 
which the evolution is followed are shown by straight lines 
in the left column. Right column: Evolution of ZBA in the 
conductance according to Eq. (|3.19|) . See the text for further 
discussion. 



a situation where the temporal variations of AsT{t) = 
Es — Et — ^coh(^) changes Tic(t) close to the above 
mentioned degeneracy point Ast = 0. In this high- 
temperature weak coupling regime the zero bias anomaly 
(ZBA) in tunnel conductance obeys the law 



-^~ln" 2 (T/T^). 

LrQ 

Substituting f3~T5|) in ([335]) . one gets 
G(t) 



Go 



(HT/Tko - vHTKo/A ST (t)y 



(3.19) 



(3.20) 



The results of the numerical analysis of the Kondo 
response to the periodical input given by Eq. 

(|3.13p as well as the adiabatically varying conductance 
([3.19P are presented in Fig. [3] It is seen from this 
analysis that the input sinusoidal potential v g cos fit 
transforms into periodic oscillations of the Kondo-related 
ZBA in G(t). The sinusoid is reproduced with a slight 
distortion when the gap Asr(t) remains positive under 
temporal perturbation (row 1 in Fig. [3]). Additional 
minima impart the " Kremlin wall" shape to the periodic 
curve G(t) when the sign of Agr(t) changes under the 
periodic perturbation (rows 2,3). The same regime with 
larger amplitude v g may result in complete suppression 
of Kondo tunneling due to the periodical triplet-singlet 
crossover (row 4). Finally, if the system remains 
completely in the singlet sector Asr(t) < near the 
crossover point, the charge perturbation results in a 
pulsed Kondo output signal (row 5). 
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It is worth noticing that this mechanism of adiabatic 
transformation of a charge signal into a Kondo response 
is close to that proposed for Kondo shuttling^ where the 
source of time-dependence are the nanoelectromechanical 
oscillations of a quantum dot with even occupation be- 
tween two leads. 

The reference Kondo temperature Tko is given by the 
equation 



T A 'o = De bo+i s ) (3.21) 

where j A = ppJ A , pp is the density of electron states in 
the leads at the Fermi level, and D is the characteristic 
scale of these states in SW regime.— This temperature 
also oscillates adiabatically in time due to the correction 
(|3.17[) to the second term in the exponent. One may 
estimate adiabatic oscillations ST^ in the lowest order 
in j'f d /(J(f + Jq) < 1- One derives from (|3~2Tj) 

Jad 

oT KO (t) ~ T K0 , T s , 2 ~ +K0—pr- ■ [3.22) 
Uo "T" JO I Jo 

These temporal oscillations only weakly distort the 
curves shown in Fig. [3] because the Kondo tempera- 
ture in these curves changes by the order of its mag- 
nitude due to oscillations changing the sign of Asr(t)- 
One may roughly estimate this effect by averaging (|3.22p 
over an oscillation period. Like in the Kondo shuttling, 15 
this averaging results in an effective enhancement of Tko 
which resembles the Debye- Waller enhancement of neu- 
tron scattering intensity because of a nonzero average 
quadratic displacement induced by phonon vibrations. 



B. Stochastic input 

We begin with the discussion of the influence of an in- 
coherent input on the spin state of quantum dot isolated 
from a metallic reservoir. As is known from the general 
theory of dynamical symmetries^ only those states from 
the total manifold are involved in its formation whose 
energies are comparable with the energy scale of the in- 
teraction which breaks the symmetry of the Hamiltonian. 
In our case this scale is determined by the Kondo tem- 
perature Tk ~ Ast ^ Apr- 



1. Fluctuations of the global constraint 

The mechanism of conversion of the stochastic com- 
ponent Sv g (t) of the input signal into a stochastic spin 
response is quite unusual. Instead of dephasing due to 
time-dependent spin flip processes,— stochastization of 
the energy spectrum of DQD results in the loss of a Curie- 
type spin response at some characteristic energy £. This 
effect is related to the time dependence of the factors 



<5stoch.A(£) (|3.8p in the Hamiltonian (|3 . 6|) . Indeed, insert- 
ing (|3.7p into (|3.6p . one may write the stochastic part of 
Hdot in the form 

ffdot* = [<Wh,s(*)P 2 - 6 stoch:E (t)M 2 }/2 (3.23) 

Unlike the adiabatic part of time dependent energies 
E\(t), this term describes spin fluctuations related to 
the dynamical symmetry of DQD. In the reduced sin- 
glet/triplet subspace the exact confinement preserved by 
the last term in (|3.6[) obeying 50(5) dynamical symme- 
try transforms into fluctuating confinement in the effec- 
tive Hamiltonian, 

H dot (t) = ^(E T S 2 + E S P 2 ) (3.24) 
- /1 (S 2 + P 2 -3) + Ws(t)P 2 

where the SO(A) symmetry is preserved only approxi- 
mately. Thus the stochastic component of the dot Hamil- 
tonian given by the correlation functions (|3.8|) appears 
explicitly in the constraint. 

It follows from (|3.24p that this component survives 
even in the asymptotic regime, T <C Ast, where the 
T/S excitations are quenched in the Kondo scattering, 
but the singlet component of the spin multiplct still in- 
fluences the constraint via its stochastic constituent. The 
effective dot Hamiltonian in this limit has the form 

H dot {t) = \e t S 2 - p(S 2 - 2) - S stoch>s (t)S 2 . (3.25) 

so that the fluctuations of the charge S/E gap may be 
transformed into the fluctuations of the spin constraint. 
This unusual situation is considered below in greater de- 
tail. 

To investigate the influence of stochastic corrections of 
jj, on the spin properties of isolated dot, we rewrite the 
Hamiltonian (|3.25[) in a fcrmionized form 

H dot (t)= ]T [E T /2 - ii{t)]ftf v . (3.26) 

u=0, ±1 

Here f v are spin fermions representing the S=l triplet,— £ 
and the time-dependent chemical potential for spin 
fermions is defined as fi(t) = /io~^stoch,s(^)- The stochas- 
tic component of p may be treated as a random potential 
in the time domain, which describes the fluctuations of 
global fermionic constraint 4 The problem of propagation 
of spin fermions in a random time-dependent potential 
may be considered by means of the "cross technique"— 
developed for the study of electron propagation in a field 
of impurities randomly distributed in real space. 

In Ref. H the short-time (white-noise) fluctuations of 
a global fermionic constraint have been considered. The 
noise correlation in this limit is delta-like 

D(t - f) = h 2 (ii(t)Kt')) - r S(t - t'), (3.27) 

Such a description presumes that the chemical potential 
suddenly " shaken" at any moment does not keep memory 



7 



about its previous value (correlation time equals zero). 
The decolierence time calculated in Born approximation 
is given by 

ft/ U ~ r 

Here we propose another realization of the stochastic 
potential, which corresponds to the situation when the 
chemical potential varies slowly in time (~ exp(— jt)). 
A very long relaxation time r r ~ I/7 with small 7 is 
assumed, so that the noise correlation is given by 

D(u) = lim ^ „ = 2tt( 2 S(lo) (3.28) 

In this limit the averaged spin propagator describes the 
ensemble of states with chemical potential fi = const 
in a given state, but this constant is random in each 
realization^ 

The problem of decoherence of the spin state in 
a stochastically perturbed DQD in this limit can be 
mapped on the so-called Keldysh mode l 26 ' 27 i 28 originally 
formulated for systems which arc ^-correlated in the mo- 
mentum space impurity scattering potential. The prob- 
lem can be solved exactly and the decoherence time is 
defined by the variance £ 2 of the Gaussian correlation 
(see below). We look for a solution of the time-dependent 
model where time is the only current coordinate in the 
system. 

The spin-fermion propagator at T = is defined as 
G&(t-f ) - (f»(t)P u (t')} R = -i{[f u (t)ti(t')} + ). (3.29) 

We sum the perturbation series for the Fourier transform 
of this Green's function 



G R (e)=g(s) 



(e) 



71=1 



(3.30) 



Here g(e) = (s + iS)^ 1 is the free spin-fermion propaga- 
tor with Et/2 — /Jo = taken as the reference energy. 
The index v is omitted, since the fluctuations of fi are 
related to the global U(l) symmetry. The noise corre- 
lation function (|3.28|) is normalized in such a way that 
corresponding vertices are dimcnsionless. The Feynman 



c) 



d) 



FIG. 4: (a) Feynman diagrams for the self energy with vertex 
corrections; (b) Bare self energy; (c) First line correction; (d) 
First vertex correction. 

diagrams for the self-energy S(e) are shown in Fig. 2] 



The key features of the Keldysh model stem from the fact 
that all self energy diagrams in a given order are equiva- 
lent due to the delta- function character of the correlation 
function D(u>) (|3.28[) . As a result, the contribution of all 
diagrams in a given order n is completely determined by 
the combinatorial coefficient A n = (2n — 1)!! giving the 
total number of diagrams corresponding to all possible 
pairwisc connections of n vertices by wavy lines. Then 
the exact analytical equation for the self energy may be 
derived^^ 



S(e)= / ^T(e,e- U ;u)G(e-u)D(u>) 



(3.31) 



Here T is the full vertex (triangle), G is the full Green 
function (thick line) and D is the noise correlation func- 
tion (wavy line) in Fig. [4jt. Evaluation of the integral 
(pPTj) with the S- functional D(w) (|3~2"8"]) gives 



E(e) =C 2 r(e,e;0)G(e) 



(3.32) 



In order to find G(e) we use the Ward identity illustrated 
by Fig. [5]) which connects the triangular vertex and the 
Green's function (GF) 



r(e,e;0) 



dG-^e) 
de 



(3.33) 



Then the Dyson equation for the spin-fermion propaga- 






FIG. 5: Diagrams for the vertex V. 



tor is transformed with the help of (|3.33p into an ordinary 
differential equation: 



dG 

~d7 



eG-l = 



(3.34) 



This equation is supplemented by the boundary 
conditio n 27 i 28 



G(e^oo) = - 



(3.35) 



The solution of (|3.34[) satisfying the boundary condition 
(|3.35[) is given by 



G R (e) 



-z 2 /2C 2 



dz 



t — z + id 



(3.36) 



Remarkably, the spin-fermion GF in this model has no 
poles, singularities or branch cuts. The solution (|3.36|) 
represents the set of spin states under a stochastically 
fluctuating chemical potential averaged with a Gaussian 
exponent characterized by the variance £ 2 . 
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Let us investigate the spin response of this "stochas- 
ticized" DQD. To calculate the spin susceptibility, it is 
convenient to make an analytic continuation of G R on the 
imaginary semi- axes of complex energies, i.e. to go over 
to the thermodynamical Matsubara Green's functions: 



1 



dz 



ic, 



(3.37) 



The spin (triplet) susceptibility at finite temperatures 
defined by the diagrams of Fig. [5] may be calculated by 
means of these functions. The spin susceptibility with 





FIG. 6: Diagrams for bare and dressed spin susceptibility. 

vertex corrections is determined as 

X(iu m ) = (3.38) 



The Ward identity (|3.33[) provides us with the exact 
equation for the vertex function 



r fl (e,e;0) = 



eG R -l 
C 2 (G«) 2 ' 



(3.39) 



giving access to the exact evaluation of the static suscep- 
tibility x(0) (right diagram on Fig. Combining (|3.38[) 
with f337| and (|3l?9"l) . we find 

rf0) ^7fe£' i! ' eV/2! " f " lh (l?) (3 ' 40) 

The asymptotic behavior of the static susceptibility 
X(0) is 



x(o) = 



Cc/T, T»C 

Gk/C, C»r 



(3.41) 



where Gc and Gk are constants. At high T the behavior 
of the dot is Curie-like with Curie constant Gc modified 
by averaging. At low T the noise dispersion £ plays the 
role of an effective temperature in the Kcldysh model 
with the corresponding constant Gk in the numerator. 

There is a great simplification in the calculations of 
the dynamic susceptibility at temperatures T > (. We 
notice that r — * 1 in this limit since transferred energy 
exceeds the dispersion of the noise spectrum. Under this 
condition one can neglect the vertex corrections, and the 
spin susceptibility is given by 



2ttC 2 



sinh 



Z2 ~ Zi 

2T 



cosh(^) cosh(^) 



+ z 2 - z 1 



Performing the analytic continuation, one gets the follow- 
ing equation for the imaginary part of the spin response 
function at real frequencies: 

Imx R (u) = 



(3.42) 



e 4 < 2 



tanh 



(-) 

\2TJ 



dq 



e 4 ( z 



cosh(2^ , 



+ 1 



cosh(^ 

Thus Imx R ~ w at small lo <^T . 

The real part of the static susceptibility is given by 

Re X R {io) = 

1 f°° , usinh(^) 



(3.43) 



du- 



-e 4 ^ 



dq- 



e 4 <^ 



1 cosh( ) + cosh( ^ ) 
where the principal part of the integrals is taken. 









U! 



FIG. 7: (Color online) Static susceptibility x(0, T). The inset 
shows a frequency dependence of Imx with a maximum at 

LO ~ 



It follows from j3~4Q]), (pPSj) and (f3T43|) and from re- 
sults of numerical calculations presented in Fig. [7] that 
the low-frequency response of a DQD in the Keldysh 
regime has nothing to do with the behavior of a free spin. 
This means that in spite of the fact that at high T a DQD 
behaves like a quantum object with spin 1, it looses at 
{oj,T} ( the generic characteristics of a localized spin 
due stochastization, hence it cannot serve as a source of 
Kondo screening at low energies. 





FIG. 8: Diagrams for spin-electron loops responsible for 
the Kondo effect. Here solid and dashed lines stand for 
spin-fermion propagators Q(ie n ) and electron propagators 
g(p,iui n ), respectively. 

A direct calculation of spin-electron loops responsible 
for Kondo screening (Fig. [8]) confirms this conclusion. 
These second order perturbation theory corrections to 
the spin-electron vertex are given by 
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Substituting Green's functions into this integral, we have and the real part of the susceptibility 



dze 



D 



tanh 



2T 



This vertex correction reflects the averaging of the stan- 
dard Kondo vertex with the Gaussian distribution of 
chemical potentials. Evaluating the integral, one obtains 
a combination of logarithmic, hypergeometric and imag- 
inary error functions 



T (2 \e -> 0)/J = poJHV2CD/0+ 



(3.44) 



:PoJ 



I 3 2 ^ 
1F1 [ -,2,— 



7rErf i 



(vac) 



where 7 = InC is the Eulcr constant. In two limiting 
cases of low and high temperatures relative to the disper- 
sion £ of noise spectrum, it leads to the following compact 
expressions 



pJ\n(D/T), T>C 
pJln(D/C), C»r 



(3.45) 



We conclude from here that the variance C 2 predeter- 
mines the energy/temperature cut-off (similarly to the 
Kondo-spin glass problem^) . This result correlates with 
the above observation that at low T <g; £ the magnetic 
excitations in stochasticized DQD loose the properties of 
spin flip processes which are essential for Kondo screen- 
ing. The characteristic decoherence time is given by: 

h/r d ~ C 

If ( -C Tk, the noise effect is seen in the behavior of the 
magnetic susceptibility as a logarithmic correction, 

X(T) = Xc (1 - ln-^max^O/Tjs-) + ...) , 



but to study the influence of 5v g (|2.6p on the Kondo 
processes at finite temperatures one should also consider 
the dephasing effects. 

We comment also the S-T response function 

X f(t - t') = -i([P a (t),Pp(t')]) -> 5 a p X p{lUn) 

which corresponds to the bare loop represented in Fig. 
HJl, where one of the two lines corresponds to the singlet 
fermionic GF while the other one represents the triplet 
fermionic GF. Since the singlet line is not affected by 
the noise, only one of the GFs in the loop suffers from 
Gaussian averaging. 

The straightforward calculations lead to the following 
answer for the imaginary part 



Imxp(uj) 



(3.46) 



1 sinh(^) 

' A ST 
v IT 



C sinh 



cxp( 



K 1 



cosh( 



2T ' 



cosh( 



Agy — u > 
2T 1 



Re X %{uj) 
1 



(3.47) 



d _ (z - A ST ) sinh {{z - A ST )/2T) ^ 



(z - A ST ) 2 



_ z 
e 2c 5 



cosh(z/2T) cosh(A ST /2T) 

The static susceptibility which mimics the Curie law 
at very large temperatures T 3> (Ast,C)i is suppressed 
exponentially ~ exp(— A$t/T) at low temperatures T <C 
(Agr, C) and has an intermediate asymptotic behaviour 
Xp(0) ~ 1/CHC/T) if A ST < T < C while XP (0) ~ 
1/Asr when C < T < A ST . The real part of the 
dynamic susceptibility taken at the resonance frequency 
w = ±Ast grows as l/£ln(£/T) when the temperature 
is lowered. The exponential suppression does not oc- 
cur since A$t is compensated by the external frequency. 
We therefore conclude that the noise may strongly affect 
the non-equilibrium clcctric-ficld-induccd Kondo trans- 
port in the regime when the singlet is a ground state of 
DQD while the access to the triplet state is facilitated by 
the applied gate voltage. 



2. Fluctuations of the scattering phase 

In accordance with the general approach to decoher- 
ence and dephasing effects^ the latter phenomena arise 
due to scattering processes at finite energy and/or tem- 
perature. These processes are described by the effective 
cotunneling Hamiltonian (|3 . 1 1 1) . The main contribution 
to dephasing is given by the term 



H, 



(stoch) 



J. 



stoch 



(t)p. 



(3.48) 



Here ^ s t oc h(^) ^ s the stochastic component of the indirect 
exchange integral calculated in Appendix A (Eq. IA3[) . 
To reveal the dephasing mechanism, one should notice 
that the parameter J^ t (t) 



w 



2V Vg(t)Mt) 



e 2 - M s Aes 



(2 



Mt) 2 (3.49) 



is in fact the modification of the effective SW exchange 
due to temporal fluctuations of the intradot exchange 
2aV, which is nothing but the gap A$t (see Eq. 13. ip . 
Unlike the similar term 6 s t oc h(t) in Eqs. (|3.7p . (|3.24[) . it 
does not influence the modulus of the vector P, but the 
components of this vector, and thereby it affects the com- 
ponents of the spin vector S via scattering processes illus- 
trated by Fig. [HI The time-dependent exchange vertex 
(|3.49p is taken in the form J so (p s (t), so that the fluctu- 
ating part in parenthesis is represented by a single mode 
¥>(*)• 

The wavy lines in the diagram stand for the correla- 
tion functions S(t — t') = (ip(t)(p(t')}, the solid line in 
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FIG. 9: Left: bare vertex in the fluctuation part of the ex- 
change Hamiltonian (|3.48[) . Right: first non-adiabatic correc- 
tion to the exchange vertex J in (13,111 1. 



the vertex correction corresponds to the bare retarded 
propagator g s (t — t') = (f s (t)f^(t'))n for pseudofermion 
excitations representing the singlet mode in the 50(5) 
multipleti2«li As was noted in I, the eventual source of 
dephasing is the gauge fluctuation field induced by non- 
adiabatic excitations in the non-diagonal operator P with 
components P v = /J/„. 

The problem of dephasing due to slow fluctuations 
S(u> — > 0) was analyzed in I, so we do not repeat here 
the corresponding calculations. The net result of this 
analysis is that the dephasing processes are relevant at 
high enough temperatures, and at T <C A.st dephasing 
is effectively quenched. 



IV. TRIPLE QUANTUM DOT IN A 
TRIANGULAR GEOMETRY. ODD ELECTRON 
OCCUPATION 

In this section we sort out the coherent and stochastic 
components of the weak probe time-dependent potential 
v g (t) (|2.6[) applied to the triangular TQD shown in Fig. 
[T}d. To make the mathematical treatment more transpar- 
ent, we consider the three-terminal geometry, so that the 
system possesses a perfect triangular symmetry C^y. It 
is convenient to enumerate both the dots and the leads 
by the numbers 1,2,3. The bias is supposed to be created 
between the leads 2,3, and the role of a passive termi- 
nal 1 is to serve as a reservoir for Kondo screening of the 
electron spin in the dot. The gate voltage v g (t) is applied 
to one of the electrodes forming the dot 1. 

Then the band Hamiltonian has the form 

Hband = E £ 3kC f jkcr C :jka . , (4.1) 
3=1,2,3 

Correspondingly, the tunnel Hamiltonian is written as 
H tun = W^2(cl k(T d jr7 + H.c.), (4.2) 

jka 

We study the excitation spectrum of the TQD in the 
charge sector Af = 1 . The spin degrees of freedom obey 
the SU(2) symmetry, and all dynamical symmetry effects 
are related to orbital degrees of freedom. We will show 
here how the influence of charge input on the orbital de- 
grees of freedom may be converted into Kondo response. 
The spectrum of the TQD was discussed in Refs. H0II22I . 



This dot is described by the Hamiltonian 

3 

Hill = e E E d Ui° + Q E ( 4 - 3 ) 

3=1 a 3 

+o' E E + v E E(4^ + H - c -)- 

(ji) ° ui) ° 

Here n jc = d] a d ]a , (jl) = (12), (23), (31), Q and Q' are 
intradot and interdot charging energies (Q ^> Q'), V is 
the interdot tunneling amplitude. In the case Af = 1 
charging terms are irrelevant and t|4. 3|) is easily diagonal- 
ized 

^ = E £ r4^r CT (4-4) 
r,cr 

Here the index F = A, E± stands for irreducible repre- 
sentations of the symmetry group of equilateral triangle. 
The basis of this representation is given by the cigenfunc- 
tions 

< CT = (4 CT + 4 CT +4 CT )/V3, (4-5) 
4 ±)CT = (dl + e ±2 ^dl + e^dl)/V3 ; 

(if = 27r/3) The ground state is the orbital singlet A, 
and the orbital degrees are quenched at low temperatures 
T< V. 

The triangular geometry provides us with a new pos- 
sibility for controlling quantum tunneling^. The tunnel 
current may be driven by means of an external magnetic 
field oriented normally to the plane of triangle, because 
the electron spectrum of electrons in the TQD is a func- 
tion of the magnetic flux $ through the triangle in such 
a geometry. As a result orbital degrees of freedom may 
be activated at finite magnetic field, and the possibil- 
ity opens to realize a charge-spin conversion mechanism 
already in the case of the TQD occupied by a single elec- 
tron (Af =1). The electron energy spectrum at finite <!> 
is 

er($) = e-2Vcosfp-|Y (4.6) 

and the values of p = 0, 27r/3, 47r/3 correspond respec- 
tively to r = A, E + , E-. At zero $ the states Ee± form 
degenerate doublet excitations. At finite $ the electron 
acquires chirality, the doublet is split, the levels evolve 
in accordance with (|4.6p and at <f>„ = (n + l/2)$ the 
ground state becomes doubly degenerate (here $0 is a 
quantum of magnetic flux, n = ±1,±2, . . .). The level 
evolution is shown in Fig. 110a . The dynamical symme- 
try SU(4) is involved in the signal transformation at <f> 
close to or in other words, at £r — £v ~ Tk for any 
pair of states T,V. Figure [TUk shows the evolution of 
the energy levels in the interval ^ <I> ^ 3$o, but the 
real periodicity of the spectrum is of course <J>o because 
at finite magnetic field all three states may be converted 
into each other by the appropriate choice of the gauge. 
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<t>0 2ct> 2<t> 



FIG. 10: (Color online) Left: Evolution of energy levels of 
a TQD in a magnetic field. Right: Zoomed level crossing 
around $ = 3$o/2. 

To make the notation symmetric we will consider below 
the sector <£> ^ $ < 2$ shown in Figure ITOTb). In 
this sector two states E± become nearly degenerate at 
$ approaching 3$o/2. When the symmetry of TQD is 
perfect, the point <& = 3<&o/2 is a point of accidental 
degeneracy, in which the symmetry of TQD is SU (4)^. 
This symmetry is broken by the potential v g (t). As a 
result level crossing transforms into anticrossing (dashed 
lines in Fig. [TU1 left panel). 

In order to calculate the contribution of 6Hdot(t) in 
the charge-spin conversion we need the matrix elements 
of the operators T[f ] (|2~T2"|) in the basis |r) of the C 3v 
point group, in particular those involved in the above 
level crossing/ anticrossing in the sector 3>o ^ $ ^ 2$o- 
The relevant matrix elements in the subspace \E±) arc 

(E ± \H<U\E±) = -2V sin ^(t), 

<£ ± |ff£K>=0 (4-7) 

and 

(E ± \h£\\E ± )=-V cos <ptf(t), 

(E ± \H^ t \E T ) = -Ve i ^j(t) (4.8) 

It follows from (|4.7|) . that the first-order term gives a 
purely adiabatic contribution ~ 4>i(t) (|2 . 14|) . Second- 
order corrections (first line in Eq. I4.8[) slightly change the 

adiabatic renormalization. The stochastic signal ~ <j>\ (t) 
arises from the off-diagonal matrix elements (second line 
in Eq. [48)) . 

A. Coherent input signal 

Following Ref. [l(J we use the irreducible representa- 
tions {A 7 E±} not only for dot eigenstates (|4.5|) but also 
for lead states 

4, CT = E( C I fcCT + 4<r + 4ka)/^ . (4-9) 

k 

4 ± „ = E(4. + e±2iv 4. + Ws ■ 

in the three-terminal geometry. 

In zero magnetic field the ground state is degenerate, 
the singly occupied TQD works as an effective spin 1/2, 



so that the effective spin Hamiltonian has the standard 
form JS • s^^4, where only fully symmetric combinations 
with r = A of lead and dot electrons are represented. Or- 
bital degrees of freedom become relevant near the cross- 
ing points, so that the effective Kondo Hamiltonian has 
the form 

Hcotun = E Jrr'SiT'Sr'T + JoTf (4.10) 

IT' 

in the representation (|4.5p . (|4.9[) . The spin op- 
erator for lead electrons is determined as Sp r , = 
Skk' c r kcr^ c r',k'cr'- Due to the orbital degeneracy 
(E + = E-), one more vector, namely the pseudospin 
vector T defined as 

T+ = ^\E+,a)(E-,a\, T~ = [T+\\ (4.11) 

r = l]T(|£ + ,a)(£ + ,aM£-^)<£-H)- 

together with its counterpart for lead electrons is involved 
in Kondo tunneling. There five vectors provide 15 gen- 
erators of the SU(4) group. The effective Kondo Hamil- 
tonian consists of 6 terms with corresponding exchange 
vertices, three of these vertices are relevant (including 
that for pseudospin interaction), and the corresponding 

( F*\ 

Kondo temperature T% for A± = exceeds the zero- 
field Kondo temperature by a factor of five (see Ref. 
fTol for details). Similarly to the case of ST degeneracy, 
deviation from the level-crossing point results in a sharp 
decrease of Tk, although in this case the peak is symmet- 
ric relative to the zero gap point (Inset in the left panel 
of Fig. ITT|) . The width of this peak may be estimated as 

1 K ■ 

In the vicinity of the level crossing point the E± or- 
bital components of the spin operators in (|4. 10[) are in- 
volved in Kondo screening, provided the difference A± = 
\E+ - £L| ~ T K . Like in the case of a DQD ([535]) . the 
Kondo temperature itself depends on the level distance, 
Tk = Tjc(A±). We work in the adiabatic regime and 
incorporate matrix elements (|4.7|) . (|4.8p in the energy 
terms e± (|4.6[) . We conclude from these equations and 
from Fig. [10] that only the non-diagonal matrix elements 
(E±\H^J t \E T ) are relevant. Mixing of two branches ~ <f>\ 
results in a time-dependent lifting of degeneracy of the 
orbital doublet at zero magnetic field, $ = 0, and at 
<f> f» 3$o/2- In agreement with the general rule for the 
Kondo effect in presence of dynamical symmctr y 10 i 20 i 35 , 
Tk (A±) is maximum for A± = 0, but unlike the case of 
DQD (Fig. [3]), the curve Tk-(A±) is symmetric around 
its maximum (see Fig. II H inset to the right panel). Now 
we may repeat the procedure proposed for the DQD in 
the previous section and calculate the time-dependent re- 
sponse of the ZBA peak in tunnel conductance to a coher- 
ent input signal (|4.8p . Two types of adiabatic temporal 
oscillations Tk (t) are presented in the left panel of Fig. 
1111 If A± oscillates between zero and some maxima due 
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FIG. 11: (Color online) Left panel: Time dependent Tk cor- 
responding to the evolution of A± in TQD. Inset: Tk as a 
function of A±. The intervals of this evolution are shown 
by straight lines. Right panel: Time dependent ZBA in the 
conductance in accordance with the evolution of Tk- 

to fluctuations of the gate voltage, Tj((t) reaches its max- 
imum at Qt = 27rn. The tunnel conductance G(t) shown 
in the right panel evolves with the same periodicity in ac- 
cordance with Eq. (|3.19p . If A± oscillates symmetrically 
around zero value, the period is halved, and G(i) reaches 
minima at Ot = im. In the general case when A± varies 
between — A a and +A&, the oscillations of G(t) are pe- 
riodic but not monochromatic. Experimentally one may 
turn from one regime to another by changing the mag- 
netic field (shifting the value of $) in the vicinity of the 
point $ = 3$ /2. 

We see that the situation with the SU(2) -> 577(4) -> 
SU(2) crossover is close to the case of the SO (4) — > 
50(5) — > U(l) crossover discussed in the previous sec- 
tion from the point of view of the conversion of the charge 
signal into Kondo response. Due to the fact that the or- 
bital degrees of freedom are involved in the formation 
of an effective exchange, the perturbation v g {t) directly 
affects Kondo tunneling by means of a time-dependent 
lowering of the point symmetry of the triangle induced 
by the gate voltage. 

B. Incoherent input signal 

The stochastization of orbital degrees of freedom in a 
TQD is induced by the term 5H^° ch \t) (l2~T5l) . lOl) . 
The relevant part is 

A (stoch) (i) = (E-\Hf° ch \E+) = -Vfilfye* (4.12) 

These stochastic inter-level transitions are responsible for 
the fluctuation induced avoided crossing E + — E^ (dashed 
lines in the right panel of Fig. [T0|) . One may write the 
corresponding part of the Hamiltonian in terms of pscu- 
dospin operator T, as 

Hf° ch = A<f och \t)T- + Af toch \t)T+ (4.13) 

This term should be added to the effective Kondo Hamil- 
tonian of TQD (|4TTU|) . 

Thus in this case the TQD stochastic fluctuations of 
gate voltage induce random pseudospin scattering, and 
the Keldysh approximation (|3.28[) may be used for this 



type of random potential. However, unlike the DQD case, 
the scattering has a vector character, so that the fermion- 
ized Hamiltonian (|4. 13|) has the form 

#do°t ch W = £?(%1<7 T + Q*(M9l ( 4 - 14 ) 

where g(t) is a random scattering potential which stems 
from (|4.13[) . g-\ and g[ are " pseudospin- fermions" for vec- 
tor T. Only the transversal component of pseudospin 
scattering is involved in the stochastic perturbation. 

Then, following the pattern of the scalar model 
p.28p . we introduce the correlation function C(t — t') = 
(g(t)g*(t')) and its Fourier transform ~ Gaussian vari- 
ance £, 

C(u>) = lim -4^—^ = 4< 2 <5M (4.15) 

The expansion of the Fourier transform of the Green's 
function for pseudospin operators 

F?(t-t>) = (g a (t)gUt')} R (4.16) 

[cf. Eq. (|3.29[) ] has the same form as in the scalar 
Keldysh model ) 26 i 28 , namely, 

oo 

F(e) = f(e) + B n {V20 2n f n+1 {e). (4.17) 

n=l 

Here /(e) = (e + iS)^ 1 is the free spin-fcrmion propa- 
gator, B n is the total number of 2n-th order diagrams. 
The indices R, a are omitted here and below for the sake 
of brevity. The main advantage of the Keldysh model, 
namely the equivalence of all diagrams corresponding to 
various combinations of noise correlation functions in the 
self energy (Fig. 2]) is still available. However, the essen- 
tial part of the diagrams in E(e) disappears due to the 
kinematic constraint T + T + = (or p^gj. = 0). This 
means that only the diagrams with the pseudospin op- 
erators ordered as . . . T + T~T + T~ . . . survive in the self 
energy (Fig. 2K). 

A similar version of the cross technique in a real 
space arises in electron systems in the domain of long- 
range Gaussian fluctuations near the charge density wave 
(CDW) instability, although the physical mechanism 
is radically different (alternating incoming and outgo- 
ing Umklapp fluctuations of CDW order parameter in 
1D2^ and 2D^^ systems)^ 

One may represent the diagrams for the vector Keldysh 
model in the following way. The vertices on the fermion 
line have two colors (say, black and white) correspond- 
ing to the first and second term in the time-dependent 
Hamiltonian (|4.14p . The black and white vertices alter- 
nate and the wavy lines connect only the vertices of dif- 
ferent colors. Therefore the perturbation series includes 
only the even order terms with equal number of black 
and white vertices. Following these rules the vertex cor- 
rection presented by the 4-th order diagram of Fig. rjji 
disappears. One of nonzero 6-th order diagrams is shown 
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in Fig. [12] (left). This is the first non- vanishing vertex 
correction to the self energy of the vector Keldysh model. 
The difference between scalar and vector Keldysh model 



-o 



A 



\ 



is in this combinatorial coefficient. In the scalar model 
the coefficient A n = (2n — 1)!!, and summation of the 
perturbation series results in Eq. ()3.36|) . In the vector 
model the coefficient B n = n\ due to above kinematic 
selection rules, and summation gives 



FIG. 12: First non-vanishing diagrams for the vertex (right) 
and the vertex correction to the fermionic self energy (left) in 
the vector Keldysh model. Black and white sites correspond 
to the two terms in the Hamiltonian (|4. 14f) ; pseudospin cor- 
relation functions (|4,16[1 are represented by dashed lines. 



J 



F(e) 



+ 00 



dxe~ x ' 2 / 2 ? f + °° dye-y 2 /^ 2 



(V2TT J_oo CV2^ 



\/ x 2 + y 2 + i5 e + \/ x 2 + y 2 + iS 



(4.18) 



r 



(see Appendix B for the derivation of this result). Wc 
see that Eq. (|4.18[) is the natural generalization of 
the one-dimensional Gaussian averaging (|3.36p char- 
acteristic for the scalar Keldysh model to the two- 
dimensional Gaussian averaging of a vector random 
field with purely transversal xy fluctuations. Only the 
modulus of the random field r — \/ x 2 + y 2 is aver- 
aged with the Raylcigh distribution function Vnie) = 
(e/C 2 ) exp(— e 2 /2£ 2 ), whereas the angular variable re- 
mains irrelevant due to the in-plane isotropy of the sys- 
tem. Like in the scalar model, the averaged pseudospin- 
fermion Green's function has no singularities. 

In order to find the Ward identity for the vertex and 
the differential equation for the Green's function gen- 
eralizing Eqs. (|3.33[) . (|3.34p wc propose for the vector 
Keldysh model the procedure which is an alternative to 
that used in31 for the scalar model. We calculate explic- 
itly the derivative dG^(e)/de using the same method of 
summation of the series in the r.h.s. of Eq. (|4.17[) as 
the one which was used for the calculation of the inte- 
gral representation (|4.18[) for the Green's function. This 
calculation (see Appendix B) gives the following result: 



(4.19) 



which is obviously the generalization of Eq. (|3 . 34[) with 
the similar boundary condition F(e — ► oo) = e . One 
may rewrite Eqs. (|3.34j) and (|4. 19)) in a unified way: 



eG{e) - 1 



C 2 G 2 (e)^-G-\e) 



eF(e)-l = ^ F 2 (£) jl| [eF -i (£)] j 



(4.20) 



One may check straightforwardly that the functions 
G(e) and F(e), given by Eqs. (j335]) and (JB2J) are indeed 
solutions of differential equations (|4.20p . Then, appealing 
to Eq. (|3.33[) . we define the vertex r(e, e, 0) for the vector 
model as 



(4.21) 



It is worth noting that the differential operator in the r. 
h. s. of Eq. (|4.2ip is nothing but div e in polar coordi- 
nates. This form reflects effective two-dimensionality of 
Gaussian averaging in the vector Keldysh model, which 
has been noticed already in Eq. (|4.18p . Equations (|4.19p 
and (|4.2ip facilitate the calculation of response functions 
of TQD. 

The stochastization of pseudospin manifests itself in 
the transformation of the response function x±(t) = 
(T~ (t)T + (0))n shown in Fig. [6l Now the solid lines cor- 
respond to pseudofermion propagators (i) and -Fj.(— t), 
and the vertex corrections are presented by the diagram 
[T2l (right panel) and the higher-order diagrams of that 
sort. In order to calculate the pseudospin susceptibility 
at finite temperature we address to the equation 

X(*w m ) = (4.22) 
)J r l(ie n )T(ie n ,ie n + iu) 



similar to Eq. (|3.38p . Here T a (ie n ) is the Matsubara- 
typc analytical continuation of the Fourier transform of 
the Green's function (|4.16p . Using the definition (|4.2ip 
of the vertex T and the second equation from (|4.20p , we 
express the vertex as 

ie n T{i<L n ) - 1 



r(e„,e n ,0) 



d 2 F 2 (ie n ) 



(4.23) 
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[cf. Eq. (|3.39|) ]. Following the procedure which led to 
Eqs. (|3.40[) and (|3.4ip in scalar model, we obtain for 
X_l(0) and its asymptotics the following equations 



and 



X±(0) = j 



x(o) 



tanh(|| 



(4.24) 



(4.25) 



Thus the pscudospin in the vector model looses its lo- 
cal characteristics in the same way as the spin in the 
scalar model. Accordingly, stochastization affects the 
Kondo tunneling. To estimate this effect we calculated 
the electron-pseudofermion loops similar to those shown 
in Fig. [8j but with solid lines standing for T a {ie n ). 
Since the weak coupling approach works only in the 
limit T, £ 3> A± , the pseudo-gap in the spectrum does 
not affect the logarithmic behaviour of the Kondo-loop. 
The stochastization-induced level repulsion in the vector 
model excludes states within the pseudogap both when 
T ^> £ and T -C £. In the latter case the infra-red 
Kondo cutoff is of the order of £. Therefore, the reso- 
nance Kondo tunneling in this case is controlled by spin 
degrees of freedom only in accordance with the SU(2) 
Kondo effect paradigm. 

We conclude from the above results, that the case 
of SU (4) symmetry supported by the interplay between 
spin and orbital degrees of freedom in a TQD differs rad- 
ically from the case of 5*0(4) symmetry involving only 
the spin variables. In the latter case the external charge 
noise results in the stochastization of spin degrees of free- 
dom, so that the DQD "looses" its spin moment at low 
enough energy/temperature. In the former case the spin 
1/2 is robust, and only the orbital (pseudospin) degrees 
of freedom are affected by the charge noise. Pseudospin 
stochastization means that the logarithmic divergences 
in the corresponding Kondo loops given by the diagrams 
similar to those of Fig. [5] are subject to a cut-off similar 
to that given by Eq. (|3.41|) . As a result, only the spin- 
electron loops determine the Kondo screening at low T. 
Thus a noise induced SU(4) — ► SU(2) crossover takes 
place in TQD. 



spin degrees of freedom are excited via the Casimir con- 
straint. In the case of SU(A) symmetry which determines 
the dynamics of TQD with single electron occupation, 
the signal v g (t) affects orbital (pseudospin) degrees of 
freedom, which are involved in the Kondo screening to- 
gether with conventional spin states. 

Both coherent and stochastic components of v g (t) are 
subject to a charge-spin transformation. The Kondo re- 
sponse to a coherent charge signal is close in its nature 
to oscillations of the Kondo temperature discussed in a 
context of Kondo shuttling!** The noise component in- 
troduces the stochastization of spin degrees of freedom. 
This stochastization may be complete in DQD, provided 
the variance of the Gaussian noise exceeds Tk- Then the 
low-energy cutoff results in the smearing and even the 
elimination of the Kondo-relatcd ZBA in tunneling con- 
ductance. In TQD with SU(A) symmetry only the charge 
(orbital) degrees of freedom are stochasticized. Strong 
enough noise may result in peculiar noise-induced quan- 
tum crossover SU(4) — > SU(2), which may be controlled 
experimentally by varying the noise level in the input 
signal. 



APPENDIX A 

The time-dependent SW transformation may be per- 
formed in the adiabatic approximation^ In our case it 
gives the following expressions for the coupling constants 
controlling T/S and T/E transitions^ 



J s (t) = 



J E (t) = 



W-w§(t) W{l-V/A ES (t))V2Mt) 



e 2 - M s (t) 



*ES 



(*) 



, (Al) 



W - uf (t) W(l - V/A ES (t))V2~Mt) 



£2 + 02 + M E (t) 



A ES (t) 



with 



(A2) 



M s {t) = M s - C s <f>l(t), M E {t) = M E + C s <fi{t) 



s VWV2d>i(t) E VWyftcfl it) 

w 2 \ l ) — a ' W 2 \ l ) — 



A 



ES 



A ES 



V. CONCLUDING REMARKS 

In this paper we have demonstrated that the charge-to 
spin conversion of a time dependent input signal applied 
to the gate attached to a side well of a complex quantum 
dot is possible through several mechanisms involving the 
dynamical symmetry of this nano-object. Such a possibil- 
ity arises when the multiplet involved in the dynamics of 
low-energy excitations includes both charge and spin de- 
grees of freedom. In the case of SO(5) symmetry charac- 
terizing the T-shaped DQD with even occupation, charge 
transfer singlet excitons are activated by v g (t), and the 



- Mt) - «f(*), 0f (*) - Mt) - 

4(t)=v g (t)/e 2 , n E {t)=v g {t)/{e 2 +Q 2 ). 

To second order in v g (t) Eqs. (|A1|) and (|A2[) lead to 
the averaged time-dependent corrections to the coupling 
coefficients 



J s (t) — J ad + J stoch 



(A3) 




e 2 - M s 
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and a similar expression for J E {t) where e 2 — Mg is re- 
placed by e 2 + Qi + M E . 



ududtfi. The angle-independent integral (|B2[) may be 
rewritten as 



APPENDIX B 

To derive the pseudospin-fermion GF for the vector 
Keldysh model, we generalize the original Keldysh sum- 
mation procedur o 26 i 28 to the perturbation series (|4.17p . 
Using the integral representation for the T-function, nl = 
J °° dzz n e~ z , we transform the series (|4.17|) into 

F(e) = /(e) J 1 + 2 £ tdt fc/(e)] 2 "] j 



(cf. Ref. \2% . Here we substituted t 2 for the variable z. 
Then changing the order of summation and integration, 
we transform F a {e) into the integral 



F(e) 



1 

1 



dx / dye-^ + y 2 y^ 2 

OO J —OO 

1 



\J x 1 + y 2 + iS e + yjx 2 + y 2 + iS 



(B3) 



which is in fact the expression in Eq. (|4.18|) . 

In order to calculate the derivative dF^/de, we start 
with the same expansion (|4. 1 7[) . The analog of Eq. (|B1[) 
for the derivative has the form 



dF(e) 
de 



»» a( f~,'£gr, w «-* 

1 - 2t 2 £ 2 f 2 (e) 



The subsequent variable change which gave Eq. (|B2|) for 
the GF gives for its derivative the following equation 



(Bl) 



Taking into account the explicit form of the free 
pseudospin-fermion propagator, fe) = (e + iS) _1 , we 
change the integration variable once more, t = u/\/2£ 
and put it into the form (|B1|) 



2£ 2 \_ e — it + i(5 e + u + id 



(B2) 



Next we introduce the "cartesian" coordinates, x 
u cos 0, y = sin0, so that u = \f x 2 + y 2 and dxdy 



dF{e) 
de 



1 J ± _J 1 

2 U 4 e 



(B4) 



where 



Jn = rfzz"cxp (-^) [/(e - - /(e + «)] 

After some manipulations, these integrals are represented 
via the GF for the vector model (IB2I): 



2e£ 2 F - 2^ 2 , J 4 = -C 4 + e 2 J 2 (B5) 



J 2 



Substituting these integrals in Eq. (|B4|) . we come after 
some algebra to the differential equation (|4. 1 9[) . 
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